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SOME  MODELS  RELATING  TO  ENCOUNTERS 
BETWEEN  SHIPS  IN  TRANSIT  AND  SUBMARINES  ON  PATROL 


by 

Gijsbertus  E.  van  DOMSELAAR 


ABSTRACT 


Models  are  developed  that  describe  encounters  between  convoys,  or  indepen¬ 
dent  ships,  and  submarines.  Two  aspects  of  such  encounters  are  considered 
i n  some  detai 1 : 


The  rate  at  which  encounters  occur  between  submarines  and 
their  targets 

The  time  that  a  submarine  remains  within  convoy  defences 
when  attempting  to  reach  an  attacking  position. 


INTRODUCTION 

Recently  in  force  employment  studies  several  high  level  campaign  models 
have  been  developed.  These  models  mainly  describe  the  interactions  between 
Soviet  submarines  and  NATO  reinforcement  and  resupply  shipping.  Values  of 
input  parameters  are  required,  which  could  either  be  determined  by  sub¬ 
modelling  or  simply  be  given.  The  choice  highly  depends  on  the  sensitivity 
of  the  results  to  these  input  parameters  and  to  the  available  information 
for  developing  submodels. 

This  memorandum  describes  submodels  for  two  input  parameters  to  higher 
models:  the  encounter  rate  between  transitting  NATO  ships  and  a  submarine 
on  patrol,  and  the  average  time  a  submarine  will  spend  inside  the  convoy 
defences  when  attempting  to  attack  a  convoy..-. 

/ 

In  section  1  a  model  for  the  encounter  rate  is  developed.  A  well  known 
encounter  model  appears  as  a  special  case.  In  section  2  a  model  is 
developed  for  the  average  time  a  submarine  needs  to  pass  a  convoy  screen. 
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1  ENCOUNTER  RATE 

1. 1  Detection  rate  model 

In  this  section,  a  model  for  the  detection  rate  between  shipping  in  transit 
(convoys  or  individual  ships)  and  a  submarine  on  patrol  is  derived.  The 
final  model  is  well  known  and  similar  results  can  be  found  in  literature 
(for  instance,  see  <1>).  The  derivation,  however,  is  given  as  an  introduc¬ 
tion  to  the  encounter  model,  for  which  model  it  will  show  its  usefulness. 

The  configuration  of  the  problem  is  illustrated  in  Fig.  1. 


The  patrol  area  of  the  submarine  has  length  L  and  width  W.  No  boundary 
effects  are  considered,  therefore  it  holds  that  d  is  much  smaller  than 
both  W  and  L  . 

The  targets  (convoys  or  individual  ships)  are  entering  on  a  side  W  and 
the  point  of  entrance  is  uniformly  distributed  over  W  .  The  targets  are 
transitting  through  the  area  on  a  straight  line  with  a  speed  u  .  The 
submarine  is  doing  a  random  search  in  the  area  with  a  speed  vg. 

The  detection  range  of  a  submarine  against  a  target  is  ^  . 


FIG.  2 
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If  <J>  is  the  angle  between  the  target  and  submarine  speed  (see  Fig.  2)  then 
the  speed  of  the  target  relative  to  the  submarine  becomes,  using  the  cosine 
rule, 


w 


J  v2 +  u2  -  2  uv  cos  <f 
■v  s  s 


(Eg.  1) 


The  submarine  will  detected  the  target  when  it  is  entering  the  detection 
circle  (radius  ^  )  around  the  target. 


FIG.  3 


During  time  interval  At  the  submarine  can  only  enter  the  detection  circle 
when  it  is  inside  the  shaded  area  of  Fig.  3.  It  is  assumed  that  the  angle 
<|»  between  u  and  v&  is  uniformly  distributed  over  the  interval  [0,  2 n]. 

Therefore  the  probability  that  the  submarine  detects  the  target  during  time 
interval  At  ,  given  the  submarine  is  within  the  angle  between  <J>  and 
(J>  d<J)  ,  becomes 


dw  d$ 
LW  2n 


At. 


(Eq-  2) 


The  rate  per  At  that  the  submarine  detects  a  target  which  is  inside  the 
patrol  area  then  becomes 


In 

1  d  J  w  d(f>  . 

2n  LW  o 


If  the  targets  are  entering  the  area  with  a  rate  of  p  per  At  ,  then  at 
any  moment  the  number  of  targets  inside  the  area  is: 


M 


L 

u 


(Eq.  3) 
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The  rate,  at  which  a  submarine  will  detect  the  targets  in  the  area  per  At, 
becomes 


2  n 

L  1  d  ,  ,,  .. 

£  u  2^  LW  ;  w  d<t>  • 
0 


Using  the  relations  cos<|>  =-  cos(n-<j>),  cos  $  =  1-  2sin2  (g) 
and  Eq.  1  the  detection  rate  becomes  ~ 


u+vs  2  71/2 


4uv_ 


£  =  »  W  T-h*  sir|2  4*  d<t> 


£  =  M  u 


^s  2  p 

W  u  n 


where  E 
second  kind. 


2  fuv^ 
u+v 


uv 


u+v  ’  2 
s 


(Eq.  4) 


2  is  the  complete  elliptic  integral  of  the 


1.1.1  Examples 

Two  special  cases,  in  which  the  targets  are  convoys,  are  derived  and  one  of 
them  is  compared  with  Eq.  4  in  a  numerical  example. 

1)  The  submarine's  search  speed  equals  the  convoy  speed, 

i.e.  u  =  vg.  Since  E(l,  ^  )  =  1,  it  follows  that  the  detection 
rate  becomes: 


e 


1 


P 


d 

W 


2)  Suppose  that  the  submarine  is  stationary  somewhere  in 
its  patrol  area,  i.e.  vg  =  0.  In  that  case  k  =  0 

and  knowing  that  E(0,  ^  ,  it  yields  the  trivial 

detection  rate 


-e2  =  M 


d 

W 


(Eq.  5) 


The  expression  for  e2  has  often  been  used  as  an  approximation  of  the 

detection  rate.  For  an  example,  the  values  of  e  (Eq.  4)  and  &2  ^ave 
been  compared.  The  following  input  values  have  been  taken: 
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u  =  15  kn 
d  =  200  n.mi 
W  =  300  n.mi 
p  =0.33  per  day. 

Figure  4  shows  the  graph  of  e  for  different  values  of  the  submarine 
search  speed  vg  and  the  graph  of  £ ^  •  The  fiQure  demonstrates  the 

increasing  detection  rate  for  increasing  submarine  speed. 


SACLANTCEN  SM-145 


1. 2  Encounter  rate  model 

After  a  detection  a  submarine  will  not  always  encounter  the  target. 
Because  of  speed  limitations  it  is  often  not  able  to  do  so.  But  also  in 
the  case  when  the  submarine  is  able  to  encounter  a  target,  it  will  not 
necessarily  try  to  intercept  the  target.  When  the  submarine  is  confronted 
with  plenty  of  targets  it  could  select  only  the  high  value  ones. 

However,  in  this  section  it  is  assumed  that  the  submarine  will  always  aim 
for  the  target  when  it  can  be  intercepted.  The  submarine  will  change  its 
course  such  that  the  relative  speed  is  towards  the  target  and  that  it  will 
change  its  search  speed  v  to  a  closing  speed  v  . 


FIG.  5 


When  the  closing  speed  of  the  submarine  is  superior  to  the  speed  of  the 
target,  the  target  will  always  be  encountered  after  a  detection  and  the 
encounter  rate  equals  the  detection  rate  of  Eq.  4: 

A  =  e  . 

When  the  submarine's  closing  speed  is  lower  than  the  target's  speed,  the 
submarine  will  only  be  able  to  encounter  the  target  when  it  is  entering  the 
detection  circle  within  the  well  known  limited  angle  of  approach: 


vc  v 

-  arcsin(  — )  <  ijj  <  arcsin(^p)  ,  (Eq.  6) 

where  iji  is  the  angle  between  u  and  w  . 

This  angle  is  illustrated  in  Fig.  5  with  the  thick  line  on  the  edge  of  the 
circle.  Within  this  angle  the  submarine  is  able  to  reach  the  centre  of  the 
target.  However  in  reality  the  submarine  does  not  need  to  reach  the  centre 
of  the  target  but  only  to  come  within  a  certain  range  of  the  centre  of  the 
target  in  order  to  fire  its  torpedoes  or  short  range  missiles.  This  range 
will  be  called  the  weapon  firing  range  r  .  From  Fig.  5  it  can  be  seen 
that  in  this  case  the  limited  angle  of  approach  becomes 

-  i(>*  <  i[i  <  4>*  , 


7 
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where 


4<*  =  arcsin  (-^)  +  arctg  ) 


(Eq.  7) 


Figure  6  shows  that  only  when  the  submarine  is  within  the  shaded  area 
during  time  interval  At  ,  will  it  be  able  to  encounter  the  target. 


Again  assuming  that  the  angle  <|>  between  u  and  v  is  uniformly  distri¬ 
buted  over  [0,2n],  the  encounter  rate  becomes 


A  =  2n  h  -f  d*  w  d<J>  .  (Eq.  8) 

o 

Obviously  the  value  of  d*  depends  on  the  value  of  the  angle  i|>  .  In 
order  to  calculate  d*  three  cases  have  to  be  separated.  In  Appendix  A 
the  encounter  rates  for  the  three  cases  are  derived.  Here  we  will  confine 
ourselves  to  the  results. 


1)  If  it  holds  that  (see  Fig.  6) 


v 

0  S  arcsi  n(  ~)  <  ^  -  ijj*  , 


then 


\ 

i 


=p  sin  4f* 


(Eq.  9) 
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If  it  holds  that  (see  Fig.  8) 

.  *  x  71 

4<  >o- 


(Eq.  11) 


where  <t>-^  and  are  solutions  of  the  quadratic  equation 


cos («li*  -  |)  = 


u-vg  cos# 


Note  that  at  the  boundaries 


arcsin(— j)  =  |  -  i|i* 


v  -  -? 


it  holds  that  A^  =  A2  anc*  ^2  =  ^3  ’  resPectl vely. 

It  has  been  said  that,  when  the  submarine  closing  speed  is  superior  to  the 
target  speed,  the  encounter  rate  equals  the  detection  rate.  But  it  can  be 
verified  that  at  the  boundary 


arcsi n(~  )  =  ^  ,  i.e.  vc  =  u, 


the  encounter  rate  A3  does  not  equal  the  detection  rate.  The  reason  for 

this  deviation  is  illustrated  in  Fig.  9,  in  which  figure  u  =  vc  =  vs  and 
the  limited  angle  of  approach  is 


di*  =  5  +  arctan  &-)  . 
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If  the  submarine  is  in  the  shaded  area  of  Fig.  9  during  time  interval  At  , 
it  will  detect  the  target.  But  it  will  not  be  able  to  encounter  the 
target,  because  it  is  entering  the  detection  circle  outside  the  limited 
angle  of  approach. 

The  deviation  at  u  =  v  will  also  appear  in  a  numerical  example  in 
Sect.  1.2.2.  c 


1.2.1  Special  case 

A  special  case  of  the  encounter  model  can  be  derived  when  the  weapon  firing 
range  is  taken  r  =  0. 

2r 

Since  arctg  (-^-)  =  0  ,  the  limiting  angle  of  approach  simply  becomes 

vc 

i)>*  =  arcsin(-^)  .  (Eq.  12) 

Now,  only  two  cases  have  to  be  separated: 

1)  If  it  holds  that 

v  v 

0  $  arcsin(-^)  -  arcsin(-jj) 

then 

v 

=  m  2  r  •  (Eq-  13) 


n 


J 
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This  expression  for  the  encounter  rate  is  well  known  and  often  used  to 
describe  the  encounter  rate  between  a  submarine  on  patrol  and  convoys  or 
ships  in  transit.  But  it  should  be  noted  that  the  expression  is  only  valid 
for  the  given  restrictions. 

2)  If  it  holds  that 


v  v 

^  -  arcsin(^p)  s  arcsin(-jj-)  <  ^ 


then 


^2  = 


u  d  !c  f  *2“*1  \ 
71  W  |_  u  \  n  2  J 


where 


u+vs  "'♦A  J2^S 

u  )E^u+vs>  2/'E[  u*vs  ’  2j) 


(Eq.  14) 


=  arccos  (  - 


=  arccos 


(v  *■ 

H  -  _£_  + 

vs  uvs 

h~- 

\  v 

V  s 


uvs  s 


+  V 


-  4 


-  J  V  2  +  V  2  -  U2‘  \ 

""  s  c  I 

An  example  of  this  special  case  is  given  in  the  next  section. 


uv  uv 
s  s 


1.2.2  Example 

In  this  section  an  application  of  the  encounter  model  Eq.  9  to  Eq.  11, 
where  r  >  0,  and  an  application  of  the  model  Eqs.  13  and  14,  where  r  =  0 
are  given.  Both  use  the  same  inputs  as  to  the  numerical  example  in  section 
1.1.1. 

In  the  first  application  it  is  assumed  that  the  configuration  of  the  convoy 
is  such  that  the  submarine  needs  to  be  14  n.mi  from  the  centre  of  the 
convoy  in  order  to  fire  its  torpedoes.  In  both  examples  it  has  been 
assumed  that  the  search  speed  of  the  submarine  equals  the  closing  speed. 
The  results  are  given  in  Fig.  10. 


2.1  Model 


When  a  submarir 
submarine  scree 
models  describi 
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submarine  needs  to  pass  a  screen  as  an  input  parameter.  In  this  section  a 
model  is  presented  for  the  average  time  a  submarine  will  spend  inside  a 
convoy  screen. 

After  a  submarine  has  detected  a  convoy  within  the  limited  angle  of 
approach,  t|i*,  it  will  change  its  course  and  speed  such  that  the  relative 
speed,  w  ,  is  aimed  at  the  centre  of  the  convoy.  Then,  if  is  the  angle 
between  v7  and  u.  The  configuration  is  given  in  Fig.  11. 


It  is  assumed  that  the  convoy  screen  is  a  circular  band  around  the  convoy, 
with  a  thickness  th  .  Figure  11  shows  that  the  time  required  to  pass  the 
screen  simply  becomes 


T 


(Eq.  15) 


However,  to  find  the  distribution  of  the  angle  iJj  is  less  simple.  In 
Appendix  C  the  structure  of  a  model  for  x  is  given,  which  takes  into 
account  the  relation  between  the  angle  iji  and  the  submarine  search  speed 


Here,  the  point  of  entrance  at  the  detection  circle  within  the  limited 
angle  of  approach  <|j*  is  assumed  to  be  independent  of  the  search  speed. 
Suppose  that  the  position  of  the  submarine  is  uniformly  distributed  on  a 
line  ahead  of  the  convoy,  perpendicular  to  the  transit  direction  of  the 
convoy  (see  Fig.  12). 


14 
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This  assumption  is  valid  when  the  submarine  search  speed  vg  is  relatively 

small  compared  to  the  convoy  speed  u  ,  but  may  also  be  valid  when  the 
submarine  has  some  a  priori  information  about  the  arrival  of  the  convoy  and 
stations  itself  in  front  of  the  convoy.  In  these  circumstances,  the  proba¬ 
bility  density  function  (p.d.f.)  of  4>  can  be  derived  as: 

9(40  =  cos  ip  ,  -  j  $  ip  £  |  .  (Eq.  16) 


Note  that  ip  is  a  different  random  variable  from  4»  in  the  previous 
sections. 

Since  the  submarine  will  only  aim  for  the  convoy  when  it  enters  the  limited 
angle  of  approach,  the  p.d.f.  of  ip  is  required,  given  the  submarine 
enters  the  detection  circle  within  the  limited  angle  of  approach.  This 
becomes: 

g(5»)  =  ^0  ,  -  4<*  4  iP  S  <|>*  ,  (Eq.  17) 

where 

«l>*  i  _  _ 

g*  =  S  -5  cos  41  dip  . 

-4,*  t 


Since  the  submarine  is  aiming  at  the  centre  of  the  convoy  it  holds  that 
«J>*  =  arcsin  , 


v  <  u 
c 


If  the  closing  speed  vc  >  0,  then  the  p.d.f.  becomes 

g(4<)  =  \  ^  cos  ip ,  -  arcsin  ip  $  arcsin^pj 


Applying  the  cosine  rule  to  Fig.  11  gives 
v  2  =  w2  +  u2  -  2  uw  cos  ip  . 


(Eq.  18) 


From  the  two  roots  of  this  equation,  the  submarine  will  take  the  one  which 
minimizes  the  time  spent  inside  the  screen.  Therefore 

w  =  u  cos  ip  +  /vc2  -  u2  sin2  ip  .  (Eq.  19) 


15 
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Then  the  average  time  a  submarine  will  spend  inside  the  screen,  using  the 
symmetry  between  ijJ  and  -  ijJ  ,  becomes 

th  -  -  - 
t  =2  J  -=-  g  (<J»)  dijf  . 

o 

Using  the  Eqs.  18  and  19  and  the  transformation  x  =  sin(ijj)  it  can  be 
found  that 


where  v  >  0  and  v  <  u  . 
c  c 

In  order  to  complete  this  section  a  general  model  should  have  been 
developed  for  the  case  in  which  the  submarine's  closing  speed  is  superior 
to  the  convoy  speed  (v^  >  u).  However,  such  a  model  cannot  be  developed 

because  the  time  the  submarine  spends  inside  the  convoy  screen  is  highly 
dependent  on  the  tactics  of  the  submarine.  Because  of  its  speed 
superiority  the  submarine  is  in  principle  able  to  approach  from  any 
direction.  For  instance  after  the  detection  it  could  move  in  front  of  the 
convoy  and  wait.  Thus,  each  type  of  tactics  requires  a  different  speci.ic 
model  to  be  developed. 


Example 


In  Fig.  13,  Eq.  20  has  been  applied  to  a  convoy  speed  u  of  15  kn  and  a 
screen  thickness  of  50  n.mi.  The  figure  denotes  the  average  time  the 
submarine  will  spend  inside  the  screen  for  different  closing  speeds  vc. 


4 
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FIG.  13 

The  example  shows  that  the  parameter  x  is  rather  insensitive  to  the 
closing  speed  of  the  submarine. 

Note  the  asymptotic  value  of  t  in  the  case  when  the  submarines  closing 
speed  approaches  the  convoy  speed.  This  theoretically  means  that,  in  the 

case  ^  ,  a  submarine  trying  to  encounter  the  convoy,  will  stay  inside 

the  screen  for  ever  (see  also  Fig.  14). 


FIG.  H 
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CONCLUSIONS 


In  this  report  models  have  been  developed  relating  to  encounters  between 
transitting  convoys  or  independent  ships  and  submarines  on  patrol.  A 
general  model  for  the  encounter  rate  is  developed.  The  model  takes  into 
account  that  a  submarine  does  not  necessarily  move  towards  the  centre  of 
the  target  to  make  an  encounter:  to  reach  weapon  firing  range  is  a 
sufficient  requirement.  The  effect  of  weapon  firing  range  on  the  encounter 
rate  is  shown  in  an  example. 

A  model  is  given  for  the  average  time  a  submarine  will  spend  inside  the 
convoy  screen  when  it  tries  to  make  an  encounter.  An  example  indicates 
that  this  parameter  is  rather  insensitive  to  the  submarine's  closing  speed. 

The  models  presented  in  this  report  can  be  regarded  as  sub-models  of 
parameters  of  "high  level"  attrition  models.  They  could  either  be  used  to 
calculate  directly  the  values  of  these  parameters  or  serve  as  a  background 
for  the  understanding  of  them  in  the  high  level  models. 
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FIG.  A .  1 


In  order  to  calculate  the  encounter  rate  \  ,  the  value  of  d*  has  to  be 
determined.  For  this  purpose  three  cases  have  to  be  separated. 

v 

1)  If  arcsin  (-^)  <  ^  ~  i|>* ,  which  means  that  the  limited  area 
of  approach  (4)*)  lies  inside  the  detection  area,  then 

d*  =  ysin(t]>*-i|>)  +  j  sin(<|»*  + 

=  d  sin  »J>*  cos  i|>  .  (Eq.  A. 4) 

(See  also  Fig.  6) 
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Then,  using  Eqs.  8  and  A.l  and  the  symmetry  of  0  in  the  intervals  (0,  n] 
and  [n,2n],  the  encounter  rate  becomes: 

K 

A1  =  2n  uW  2d  sin^*  (u_Vs  C0S  ^  d<t>  ’ 

o 


^  =  P  g  sin  tj)* 


(Eq.  A. 5) 


2) 


If  arcsin  ^  -  iJj*  and  tjj*  <  (see  Fig.  7)  , 


then 


d*  = 


d  sin  cos  <J<  ,  if  0  ^  tji  <  ^ -  41* 

~  sin(<J/*-<j>)  +|  ,  if  arcsin 


Using  the  relations  (A.l)  and  (A. 2),  it  can  be  found  that 


u-v  cos 


d*  = 


d  si  nil)* 


,  if  0  $  <  $1  and  <j>2  <  4  tt 


u-v  cos  <f>  v 

|[l+simj)* - — - ^  cos  4>*  si n<t>] , 


if  4^  0  $  <t>2  , 


(Eq.  A. 6) 


where  and  (j>2  are  solutions  of  the  quadratic  equation 


cos  (^  = 


u-vs  cos 
w 


The  solutions  and  <t>2  are  given  in  Appendix  B.l. 

The  following  relations  hold: 

cos  <(>2  -  cos  ^  =  -  ctg  (|  )  I  sin  <j>2  -  sin  <(^1  , 


(Eq.  A. 7) 
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/  wdd>  =  2(u+v  ) 
*1 


Relation  (A. 7)  is  proved  in  Appendix  B.2. 


U  2^"s  *'*l\ 

1 U+vs  ’  2  ) 

E\u+vs  ’  2  )\ 

(Eq.  A. 8) 


With  the  aid  of  Eqs.  A. 6  to  A. 8  the  encounter  rate  of  Eq.  8  can  be  written 
as  \2  >  presented  in  Eq.  10. 


3)  If  «|»*  >  £ 


then 


d*  = 


,  if  0  ^  «1>  <  t|>* 


(  |  sin(i|>*-4<)+  |  ,  if  t|>*-  |  $  i|/  <  arcsin 


Using  the  relations  (A.l)  and  (A. 2),  it  can  be  found  that 
/  d  ,  if  0$d><d>^  and  <t>2<<Kn 


d*  = 


|[l+sin4i* 


u-v  cos  d)  v 


(Eq.  A. 9) 


^cosijj*  Sind>] ,  if  d>-j*<K<t>2  , 


where  d>^  and  <t>2  are  solutions  of  the  quadratic  equation 

U-V  COS  d> 

cos(**-  5)=  _L_ -  . 

See  Appendix  B.l  for  the  solutions  d^  and  d^- 

Applying  Eqs.  A.  7  to  A.  9  to  Eq.  8,  the  encounter  rate  K3  can  be  written 
as  presented  in  Eq.  11. 
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APPENDIX  B.l 

THE  SOLUTIONS  0i  AND  02 


The  equation 

u-v  cos  0 

s 

COS  0  =  - -  , 

^  w 

where  w  =  yu2  +  vg2  -  2uvg  cos  0  and  0  has 
written  as 

vg2  cos2  0  +2  uvs(cos20-l  )cos  0  +  u2(l-cos  2  0)  -  v^: 

This  quadratic  equation  has  the  roots 

0!  =  arccos  /—  sin2  0  +  c-— ^  f~v  2  -  u  2  sin2  0 

\vs  s  s 

02  =  arccos  (  —  sin20  -  c°^s  ^  vg2  -  u 2  sin2  0 

\  Vs  vs  b 

Note  that  the  condition  v  2  -  u  2  sin2  0^0 
Eq.  A. 3).  5 


a  given  value,  can  be 

2  cos2  0=0. 

)■ 

)• 

is  always  valid  (see 
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APPENDIX  B. 2 

THE  PROOF  OF  RELATION  (Eg.  A. 7) 

Applying  the  sine  rule  to  Fig.  B.2.1  gives  the  relation 
u  vs 

sin  "  sTm r  ’  Where  *  ”  ♦»  * 

u 

FIG.  B.2.1 


this  relation  can  be  written  as 


si  nil)  =  sin  41  cos  q>  +  sin  $  cos  q>  . 


Substituting  the  solutions  <>  and  <(>  of  Appendix  B.l  into  this  equation 
gives 


sin  <J>j  =  sin  q*  cos  qj  -  sinq>  J sin  2  qj  , 


sin  ())2  =  sin  qj  cos 


q>  +  77-sin  i|f  U 


-  sin  i  q>  . 


Using  the  relation 


cos  q>2  -  cos  q> j  -  -  2  —  cos  q>  v  ^ —  sin z  q<  , 
it  follows  that 


cos  <t>2  -  cos4>t  =  -  ctgq>  (  sin  q»2  -  sin  0j). 


Equation  A.  7  is  found  by  taking  q>  =  ^  -  qi*. 
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APPENDIX  C 

THE  STRUCTURE  OF  A  MODEL  FOR  x 


In  this  Appendix  the  structure  of  a  model  for  parameter  x  is  derived 
which  is  more  general  than  the  model  in  Sect.  2.  At  the  moment  a  submarine 
enters  the  detection  circle,  the  convoy  is  sailing  with  a  relative  speed  w 
towards  the  stationary  submarine.  This  relative  speed  is  determined  by  the 
speeds  of  the  submarine  (v  )  and  convoy  (u). 


In  Fig.  C.l  it  can  be  seen  that  the  point  of  entrance  at  the  detection 
circle  is  determined  by  the  angle  p  .  Similar  to  the  derivation  of  the 
p.d.f.  in  Eq.  17  it  can  be  found  that  the  p.d.f.  of  the  random  variable  p 
becomes 

9(P)  =  ,  0  <  p  <  7t  .  (Eq.  C.l) 


When  the  submarine  enters  the  detection  circle  within  the  limited  angle  of 
approach,  4<*  ,  it  changes  its  course  and  speed  (from  vs  to  v  )  such 
that,  similar  to  Eq.  19,  the  relative  speed  w  becomes:  J  c 

w  =  u  sin  p  +  , /  vc2  -  u  2  cos  2  p  .  (Eq.  C.2) 

In  order  to  calculate  the  average  time  a  submarine  spends  inside  the 
screen,  two  cases  have  to  be  separated.  The  separation  is  similar  to  the 
separation  made  for  the  special  case  in  Sect.  1.2.1. 

v  v 

1)  arcsin(j^)<  ^  ,  where  41*  =arcsin  (-jj  ) 
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This  means  that  the  maximum  amplitude  of  4»  less  than  ^  -  0*  (see 
also  Fig.  C. 1). 

In  this  case  the  p.d.f.  of  the  random  variable  0  ,  given 

submarine  enters  the  detection  circle  within  the  limited  angle  of 
becomes 


9,(0)=  ,  \  '  **  ’ 


where 


Then 


n/2  -ip+ip3* 


sin  0  _  sin(0+0*)  -sin(0-0*) 

2  up  ~  2 


n/2  -0-0* 
n  n/2  -0+  4** 


T  =  2 


4^  g,(P)  f(4>)  dp  d <>  , 

w  1 


o  7t/ 2 -ip  -0* 


where  f ( <t>)  is  the  p.d.f.  of  <t>  ,  which  is  the  angle  between 

As  before  it  could  be  assumed  that  0  is  uniformly  distributed,  i.t 


f(0)  =  ^  ,  0  $  0  $  2n  . 


The  relations  between  0  and  0  are  given  by 

u-v  cos  0  v  sin  0 

cos0  =  -  and  sin  0  =  - , 

x  w  w 

where  w  is  as  in  Eq.  1. 

v 

2)  arcsin  (-jj  )  §  5  -  0*  . 

Similarly  as  in  Sect.  1.2  the  integrals  are  divided  by  the  cases 
0  S  ^  -  0*  and  0  =  ^  “  0* . 

In  case  0  $  ~  -  0*  Eq.  C.4  still  holds. 


that  the 
approach, 

(Eq.  C. 3) 

(Eq.  C.4) 

(Eq.  C. 5) 

u  and  v&. 
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In  case  *i>  ^  ^  *  *1»*  ,  the  p.d.f.  of  p,  given  that  the  submarine  will 
encounter  the  convoy,  becomes 


92  (P) 


_  sin  p 

2g2* 


0CP$5  - 


where 


7t/2— tp+ip3*' 


0 


1  -  sin(4i->i)*) 
2 


The  average  time  a  submarine  will  spend  inside  the  screen,  given  the  sub¬ 
marine  will  encounter  the  convoy,  becomes 


n/ . . 

1  =  2/  f  iS-  g  (p)  f(<t>)  dp  d<>  + 

o  rr/2-rp— 


7t  n/2-4f+i|>* 

2 ;  ;  g,  (P)  f(«t>)  dp  dci.  + 

7t/ 2— ip  ~l|i* 


<t>2  «/ 2-«|»+«|»*  th 

2  ;  /  Sg-  g2  (P)  f  (4»)  dp  d$  , 

4>,  o 


(Eq.  C. 6) 


where  f(<j>)  is  as  defined  after  Eq.  C.5,  and  <|>  and  <j>2  are  the  solutions 
of  the  equation 


n  u-v  cos  $ 

cos  (2"  >)»*)=  - - -  ,  (see  Appendix  B.  1). 
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